arXiv:1501.05566v2 [hep-ex] 19 Nov 2015 


EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2015-001 

2015/11/23 


CMS-EXO-12-057 


Search for heavy Majorana neutrinos in + jets events 
in proton-proton collisions at y/s = 8 TeV 


The CMS Collaboration^ 


Abstract 


A search is performed for heavy Majorana neutrinos (N) using an event signature 
defined by two muons of the same charge and two jets (/A^Ajj). The data corre¬ 
spond to an integrated luminosity of 19.7 fb 1 of proton-proton collisions at a center- 
of-mass energy of 8 TeV, collected with the CMS detector at the CERN LHC. No ex¬ 
cess of events is observed beyond the expected standard model background and up¬ 
per limits are set on |V;(n| 2 as a function of Majorana neutrino mass m n for masses 
in the range of 40-500 GeV, where is the mixing element of the heavy neutrino 
with the standard model muon neutrino. The limits obtained are |Gn| 2 < 0.00470 
for m n = 90 GeV, IGnI 2 < 0.0123 for = 200GeV, and |Gn| 2 < 0.583 for 
»zn = 500 GeV. These results extend considerably the regions excluded by previous 
direct searches. 


Submitted to Physics Letters B 


© 2015 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license 


*See Appendix[A]for the list of collaboration members 








1 


1 Introduction 


The non-zero masses of neutrinos and mixing between flavors have been well established by 
neutrino oscillation experiments and provide evidence for physics beyond the standard model 
(SM) [III. The most compelling way of explaining the smallness of the neutrino masses is the 
"seesaw" mechanism, which can be realized in several different schemes mm. In the simplest 
model, the SM neutrino mass is given by m v ~ y 2 v 2 / rn^, where y v is a Yukawa coupling of 
v to the Higgs field, v is the Higgs vacuum expectation value in the SM, and is the mass 
of a new heavy neutrino state (N). In this scheme, N is a Majorana particle (which is its own 
antiparticle), so processes that violate lepton number conservation by two units are possible. 

In this paper we describe a search for heavy Majorana neutrinos using a phenomenological 
approach Iil2ljl9l . We follow the studies in Refs. 11711191 and consider a heavy neutrino that 
mixes with the SM muon neutrino, with m n and as free parameters of the model. Here 
is a mixing element describing the mixing between the heavy Majorana neutrino and the 
SM muon neutrino. 


Early direct searches for heavy Majorana neutrinos based on this model were reported by the 
L3 lEOll and DELPHI ED experiments at LEP. They searched for Z —> iy N decays, where ly is 
any SM neutrino (£ = e, y, or t), from which limits on | V„n | 2 as a function of m N can be derived 
for Majorana neutrino masses up to approximately 90 GeV. More recently, the CMS experiment 
at the CERN LHC extended the search region up to about 200GeV E2j . Several experiments 
have obtained limits in the low mass region (m^ < 5 GeV), including LHCb E3ll at the LHC. 
The searches by L3, DELPHI, and LHCb allow for a finite heavy-neutrino lifetime such that it 
decays with a vertex displaced from the interaction point, while in the search reported here it 
is assumed that the N decays with no significant displacement of the vertex. Precision elec- 
troweak measurements can be used to constrain the mixing elements, resulting in indirect 90% 
confidence level limits of V ; ,n | 2 < 0.0032 independent of heavy neutrino mass E4l . Other 
models with heavy neutrinos have also been examined. ATLAS and CMS at the LHC have 
reported limits on heavy Majorana neutrino production |25726] in the context of the Left-Right 
Symmetric Model. ATLAS have also set limits based on an effective Lagrangian approach E6l . 

We report on an updated search for the production of a heavy Majorana neutrino in proton- 
proton (pp) collisions at a center-of-mass energy of = 8 TeV at the LHC with a data set 
corresponding to an integrated luminosity of 19.7 fb 1 collected with the CMS detector. We 
assume that the heavy Majorana neutrino is produced by s-channel production of a W boson, 
which decays via W + —> Ny + . The N is assumed to be a Majorana particle, so it can decay 
via N — > W y 1 with W —> qq', resulting in a y 1 y 1 qq 7 final state. The principal Leynman 
diagram for this process is shown in Lig. [T] The charge-conjugate decay chain also contributes 
and results in a y _ y“qq / final state. The lowest order cross section for pp —> (W 21 )* —> Ny 2 at 
a center of mass energy i/s is given by 


u(s) = 



d y E \f* (*' Q 2 )/? (y> Q 2 )l 


W 


( 1 ) 


where fq are the parton distribution functions for quark q at Q 2 = s = xys, and x and y are the 
fractions of the proton momentum carried by the interacting quarks. The parton subprocess 
cross section a(s) is given by |27l : 


n(s) = 


TCCi 
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where aw is the weak coupling constant and niw is the W-boson mass. 


( 2 ) 
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2 Detector, signal simulation, and data selection 



Figure 1: The lowest order Feynman diagram for production of a heavy Majorana neutrino N. 
The charge-conjugate diagram also contributes and results in a i~i~ qq' final state. 

We search for events with two isolated muons of the same sign of electric charge and at least 
two accompanying jets. The backgrounds to such dimuon final states originate from SM pro¬ 
cesses that contain isolated prompt same-sign (SS) dimuons in the final state (e.g. WZ produc¬ 
tion), and from processes such as multijet production in which muons from b-quark decays, or 
from jets, are misidentified as isolated prompt muons. The misidentified muon background is 
more significant for low masses (m^ < 90 GeV), while at higher masses the SM prompt dimuon 
background becomes more important. 


2 Detector, signal simulation, and data selection 

The CMS detector is described in detail in Ref. 1281 . Its central feature is a superconducting 
solenoid, which provides a magnetic field of 3.8 T along the direction of the counterclockwise 
rotating beam, taken as the z axis of the detector coordinate system, with the center of the 
detector defined to be at z = 0. The azimuthal angle (p is measured in the plane perpendic¬ 
ular to the z axis, while the polar angle 9 is measured with respect to this axis. Within the 
superconducting solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal 
electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of 
a barrel and two endcap sections. Muons are measured in gas-ionization detectors embedded 
in the steel flux-return yoke outside the solenoid. Extensive forward calorimetry complements 
the coverage provided by the barrel and endcap detectors. A two-level trigger system selects 
the most interesting events for analysis. 

Muons are measured in the pseudorapidity range |//| < 2.4, defined as r] = - In [tan ((9/2)], 
with detection planes made using three technologies: drift tubes, cathode strip chambers, and 
resistive plate chambers. Charged-particle trajectories are measured in a silicon pixel and strip 
tracker covering 0 < (p < 2n in azimuth and )/| < 2.5. Matching muons to tracks measured in 
the silicon tracker results in a relative transverse momentum resolution for muons with 20 < 
pj < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the endcaps. The pj resolution in 
the barrel is better than 10% for muons with pj up to 1 TeV |29l . 

The heavy Majorana neutrino production and decay process are simulated using the leading- 
order (LO) event generator described in Ref. fl8l and implemented in ALPGEN v2.14 11301 . 
The LO cross section at i/s = 8TeV for pp —* N/G —> p fCqq' with V,,\jj 2 = 1 has a 
value of 1515 pb for = 40 GeV, dropping to 3.56 pb for m n = 100 GeV, and to 2.15 pb 

for mjsj = 500 GeV fl8l . The cross section is proportional to j V„n | 2 - We scale the leading order 
cross section by a factor of 1.34 to account for higher-order corrections, based on the next- 
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to-next-to-leading-order calculation in Refs. [31, 32 j for s-channel W' production that has the 
same production kinematics as the signal. We use the CTEQ6M parton distribution functions 
(PDFs) Il33l . Parton showering and hadronization are simulated using PYTHIA v6.4.22 J34j. 
The Monte Carlo (MC) event generator is interfaced with CMS software, where GEANT4 1351 
detector simulation, digitization of simulated electronic signals, and event reconstruction are 
performed. To ensure correct simulation of the number of additional interactions per bunch 
crossing (pileup), Monte Carlo simulated events are mixed with multiple minimum bias events 
with weights chosen using the distribution of the number of reconstructed pp interaction ver¬ 
tices observed in data. 

A dimuon trigger is used to select the signal sample. This trigger requires the presence of 
one muon reconstructed from matching tracks in the muon system and tracker with transverse 
momentum (pj) above 17GeV and a second muon with pj > 8GeV formed from a track in 
the tracker matched to hits in one of the muon detectors. The trigger efficiency for the signal 
sample is measured using Z —> p + p events selected in data, and is found to be in the range 
85-95%, varying with pj and //. 

Additional selections are performed to ensure the presence of well-identified muons and jets. 
Events are first required to have a well-reconstructed pp interaction vertex (primary vertex) 
identified as the reconstructed vertex with the largest value of £ Pj for its associated charged 
tracks reconstructed in the tracking detectors C36l . 

Leptons, jets, and missing transverse energy (Et) in the event are reconstructed using the stan¬ 
dard CMS particle-flow techniques Il37l(38 |. The missing transverse momentum vector is de¬ 
fined as the projection on the plane perpendicular to the beams of the negative vector sum 
of the momenta of all reconstructed particles in an event. Its magnitude is referred to as Et- 
Jets are formed from clusters based on the anti-kj algorithm H39l . with a distance parameter 
of 0.5, and are required to be within the pseudorapidity range |?/| < 2.5. Muon candidates 
are required to have \t]\ < 2.4 and to be consistent with originating from the primary vertex. 
Muon candidates are reconstructed by matching tracks in the silicon tracker to hits in the outer 
muon system, and are also required to satisfy specific track quality and calorimeter deposition 
requirements [29[ . Muon candidates must be isolated from other activity in the event, which is 
insured by requiring their relative isolation parameter (I re i) to be less than 0.05. Here I re \ is de¬ 
fined as the scalar sum of transverse energy present within A R = \/(At]) 2 + (A cp) 2 < 0.3 of the 
candidate's direction, excluding the candidate itself, divided by the muon candidate's trans¬ 
verse momentum. The muon selection criteria are the same as those used in Ref. If40l except for 
the more stringent requirement on / re | used here. 

Events are required to contain two same-sign muons, one with pj > 20 GeV and the other with 
pj > 15 GeV, and at least two jets with pj > 20 GeV. Events with a third muon are rejected to 
suppress diboson events such as WZ. To reduce backgrounds from top-quark decays, events in 
which at least one jet is identified as originating from a b quark are rejected, where the medium 
working point of the combined secondary vertex tagger HIM has been used. 

Further selection requirements are made based on two Majorana neutrino mass search regions. 
In the low-mass search region (« 1 n < 80 GeV) the W-boson propagator is on-shell and the 
final state system of p ± p ± op[ should have an invariant mass close to In this region the 
following selections are imposed: missing transverse energy <30 GeV; p^p^p invariant mass 
<200GeV, where the two jets chosen are those which result in m(p ± p ± ) j) closest to %; and 
dijet invariant mass of these two jets m(jj) < 120 GeV. In the low-mass region the transverse 
momenta of the muons and jets are relatively low and the overall efficiency is considerably 
smaller than at higher masses. Therefore, the invariant mass requirements in this region are 
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3 Background estimation 


chosen to be relatively loose. 

In the high-mass search region (% > 80 GeV) the W-boson propagator is off-shell but the W 
boson from the N decay is on-shell, so the W —> qqf decay should result in a dijet invariant mass 
close to %. As the decay particles are more energetic in this region, the following selections 
are imposed: missing transverse energy less than 35 GeV (relaxed to increase signal efficiency); 
dijet invariant mass 50 < m(jj) < 110 GeV, where the two jets with the invariant mass closest 
to ni\\! are chosen. In both the low-mass and high-mass search regions, the upper cutoff on the 
missing transverse energy suppresses SM background processes in which a W boson decays 
leptonically (W —» }iv), including W + jet and tt production. 

The low-mass selection gives the best sensitivity for masses below 90 GeV, while for > 
90 GeV the high-mass selection gives the best sensitivity. Therefore, we use the low-mass selec¬ 
tion for 40 < fflfj < 90 GeV and the high-mass selection for rn ^ > 90 GeV. For masses less than 
40 GeV the overall acceptance for the Majorana neutrino signal is less than 1% and we do not 
search in this region. 

After applying all the selection criteria above, a final selection is applied based on optimizing 
the signal significance using a figure of merit f42l defined by e$/ ( a/2 + SB) with the num¬ 
ber of standard deviations a = 2 and where £s is the signal selection efficiency and SB is the 
uncertainty in the estimated background. 

At each Majorana neutrino mass point, the optimization is performed by varying the lower 
bound of three selections: the transverse momentum pj\ of the highest pj muon (or "leading" 
muon), the transverse momentum px 2 of the second muon (or "trailing" muon), and the in¬ 
variant mass The two jets used in this optimization are those selected as described 

above for the low- and high-mass regions. 

The overall signal acceptance includes trigger efficiency, geometrical acceptance, and efficien¬ 
cies of all selection criteria. The overall acceptance for heavy Majorana neutrino events ranges 
between 0.69% for = 40 GeV to 12% for = 500 GeV. The lower acceptance at low 
is due to the smaller average pj of the jets and muons in these events. The acceptance dips 
when ton approaches % because the muon pj becomes small. Above the W-boson mass the 
acceptance increases again up to 400-500 GeV at which point the boosted decay products of the 
Majorana neutrino begin to overlap. 


3 Background estimation 

There are three potential sources of same-sign dimuon backgrounds: SM sources such as WZ 
production, events resulting from misidentified muons, and opposite-sign dimuon events (e.g. 
from Z —> p 1 p~) in which the charge of one of the muons is mismeasured. The latter source 
is found to be negligible for muons with p-\ in the range of interest for this analysis, based on 
MC studies and studies with cosmic ray muons. Estimation of the remaining two sources of 
background is discussed below. 

An important background source is the irreducible background from SM production of two 
genuine isolated muons of the same sign, which can originate from sources such as WZ and 
ZZ diboson production, double W-strahlung W ± W ± qq processes, WH production, ttW pro¬ 
duction, double parton scattering (two qq' 4 W), and triboson production. These processes 
have relatively small cross sections, and are consequently estimated using MC simulations. We 
use PYTHIA to simulate ZZ and WZ production and MadGraph v5. 1.3.30 Ii43l for the remain¬ 
ing processes. 
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The second significant background source originates from events containing objects misiden- 
tified as prompt muons. These "prompt muons" originate from b-quark decays or light-quark 
or gluon jets. Examples of this background include: multijet production in which two jets are 
misidentified as muons; W(—» pv) + jets events in which one of the jets is misidentified as 
a muon; and tt decays in which one of the top-quark decays yields a prompt isolated muon 
(t —> Wb —> }iv fl b), and the other muon of same charge arises from a b-quark decay or a jet 
misidentified as an isolated prompt muon. These backgrounds are estimated using control 
samples from collision data as described below. 

To estimate the misidentified muon background we use the method described in Ref. (22|. An 
independent data sample enriched in multijet events is used to calculate the probability for a 
jet that passes minimal muon selection requirements ("loose muons") to also pass the more 
stringent requirements used to define muons selected in the heavy Majorana neutrino signal 
selection ("tight muons"). This probability is binned in pj and rj and is used as a weight in 
the calculation of the background in events that pass all the signal selections except that one or 
both muons fail the tight criteria, but pass the loose ones. 

The sample enriched in multijet production is selected by requiring a loose muon and a jet, 
resulting in events that are mostly dijet events with one jet containing a muon. Only one muon 
is allowed and upper cutoffs on missing transverse energy (f% < 20 GeV) and the transverse 
mass (nij < 25 GeV) are applied, where m T is calculated using the muon pj and Et- These 
requirements suppress contamination from W and Z boson decays. We also require the loose 
muon and jet to be separated in azimuth by A <p > 2.5. This away-side jet is used as a tag and 
the loose muon is a probe used to determine the misidentification probability. The transverse 
energy of the tag jet is an essential ingredient to calibrate the characteristics of the probe (loose 
muon). 

Loose muons are defined by relaxing the identification requirements (used to select signal 
events) as follows: the isolation requirement is relaxed from / re | < 0.05 to / re | < 0.4; the trans¬ 
verse impact parameter of the muon track is relaxed from <0.05 mm to <2 mm; the chi squared 
per degree of freedom of the muon track fit is relaxed from 10 to 50. 

The overall systematic uncertainty in the misidentified muon background is determined from 
the variation of the background estimate with respect to the isolation requirement for the loose 
muons and the pj requirement for the tagging jet. Increasing and decreasing the pj require¬ 
ment for the tag jet changes the pj spectrum of the recoiling muon in the event and is found 
to have the largest impact on the background level. The misidentified muon rate is compara¬ 
tively stable with respect to variations in loose muon isolation requirements. As a result, the 
28% overall systematic uncertainty in the misidentified muon background estimate is domi¬ 
nated by the px requirement on the tag jet. 

We evaluate the method used to estimate the background from misidentified muons by check¬ 
ing the procedure using MC simulated event samples in which the true origin of the muons, 
either from W or Z boson decays or from a quark decay, is known. The misidentification prob¬ 
abilities are obtained from multijet events and are used to estimate the misidentified muon 
backgrounds in tt, W + jets, and independent multijet events by applying the background es¬ 
timation method described above. The predicted backgrounds agree with the expectations 
within the systematic uncertainties. 

In addition, to test the validity of the background estimation method, we define two signal- 
free control regions in data. We apply the background estimation method in these regions and 
compare the result with the observed yields. The control regions in the two mass ranges are 
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4 Systematic uncertainties 


defined as follows. In both the low-mass range (40 < mn < 90 GeV) and the high-mass range 
(m N > 90 GeV) the control region is the same as the signal selection, without the final optimized 
selections but with either missing transverse energy greater than 50 GeV or one or more jets that 
are tagged as originating from a b quark. In the low-mass control region we predict a total back¬ 
ground of 51.4 ± 1.9 (stat) ± 8.3 (syst) events compared to an observed yield of 45 events, while 
in the high mass control region we predict a total background of 87.6 ± 2.5 (stat) ± 12.3 (syst) 
compared to the observed yield of 81 events. The misidentified muon background accounts 
for about 2/3 of the total background in both regions. The remaining backgrounds originate 
mainly from diboson production (approximately 25% of the total background) and Higgs bo¬ 
son production (approximately 4%). In both regions the predictions are in agreement with 
the observations and well within the systematic uncertainty, which is dominated by the 28% 
uncertainty in the misidentified muon background. The observed distributions of all relevant 
observables also agree with the predictions, within uncertainties. 


4 Systematic uncertainties 

The main sources of systematic uncertainties are associated with the background estimates. As 
described above, the overall systematic uncertainty in the misidentified muon background is 
28%. The systematic uncertainties in the normalizations of irreducible SM backgrounds are 
12% for WZ and 9% for ZZ f44l . For the other processes the uncertainty is 25%, determined by 
varying the renormalization and factorization scales from the nominal value of Q 2 to 4Q 2 and 
Q 2 /4, and following the PDF4LHC recommendations [45ll46| to estimate the uncertainty due to 
the choice of PDFs. After combining all the SM backgrounds with their respective uncertainties, 
the overall systematic uncertainty in the SM irreducible background estimate, including those 
discussed below, is 19%. 

Sources of systematic uncertainty associated with the estimates of the heavy Majorana neutrino 
signal and of the SM irreducible background are summarized in Table [l] To evaluate the un- 

Table 1: Summary of systematic uncertainties in the estimation of the heavy Majorana signal, 
SM background (SM Bkgd.), and misidentified muon background (Misid. Bkgd.) for low-mass 
selection (first number) and high mass selection (second number in square brackets). 


Source 

Signal (%) 

SM Bkgd. (%) 

Misid. Bkgd. (%) 

Background estimate from data 

— 

— 

28 [28] 

SM cross section 

— 

9-25 [9-25] 

— 

Jet energy scale 

5-12 [1-3] 

3 [5] 

— 

Jet energy resolution 

6-12 [2-3] 

7 [10] 

— 

Event pileup 

6-7 [3-6] 

7 [4] 

— 

Unclustered energy 

2 [1] 

3 [1] 

— 

Integrated luminosity 

2.6 [2.6] 

2.6 [2.6] 

— 

Muon trigger and selection 

2 [2] 

2 [2] 

— 

b-tagging 

1 [1] 

1.5 [1] 

— 

PDF 

3.5 [3.5] 

— 

— 

Renormalization/ Factorization scales 

8-10 [1-6] 

— 

— 

Signal MC statistics 

3-9 [1-4] 

— 

— 


certainty due to imperfect knowledge of the integrated luminosity H47l . jet energy scale [j48|. 
jet energy resolution Il48l , b-tagging IfiTI , muon trigger and selection efficiency, and the cross 
section for minimum bias production used in the pileup reweighting procedure in simulation, 
the input value of each parameter is changed by ±1 standard deviation from its central value. 
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Energy not clustered in the detector affects the overall missing transverse energy scale resulting 
in an uncertainty in the event yield due to the upper cutoff on Kt- Additional uncertainties in 
the heavy Majorana neutrino signal estimate arise from the choice of PDFs and renormalization 
and factorization scales used in the ALPGEN MC event generator. 


5 Results and discussion 

The selections determined by the optimization for each Majorana neutrino mass point are 
shown in Table [2] together with the overall signal acceptance. Figure [2] shows the transverse 

Table 2: Minimum thresholds on discriminating variables determined by the optimization and 
overall signal acceptance for each Majorana neutrino mass point (m^). As discussed in the text, 
the search was divided into low-mass (m^ < 80 GeV) and high-mass (m\j > 80 GeV) regions, 
and different selection criteria were used in the two regions. 


m N 

(GeV) 

m^^jj) 

(GeV) 

Pti 

(GeV) 

PT2 

(GeV) 

Acceptance 

(%) 

40 

80 

20 

15 

0.69 

50 

80 

20 

15 

0.80 

60 

80 

20 

15 

0.64 

70 

80 

20 

15 

0.26 

80 

80 

20 

15 

1.2 

90 

110 

20 

15 

1.2 

100 

120 

20 

15 

4.7 

125 

140 

25 

20 

11 

150 

160 

35 

25 

13 

175 

200 

45 

30 

15 

200 

220 

50 

35 

16 

250 

270 

75 

35 

17 

300 

290 

100 

45 

15 

350 

290 

100 

45 

16 

400 

290 

100 

45 

15 

500 

290 

100 

45 

12 


momentum distributions of the two muons, pxi and pj 2 , and the invariant mass m{ji ± y ± ])) for 
the low-mass region after all selections are applied except for the final optimization require¬ 
ments. The corresponding distributions for the high-mass region are shown in Fig. [3] 

After applying all the final optimized selections we obtain the results shown in Table [3] The 
expected signal depends on m n and | | 2 . For m N = 50GeV and |V ? ,n| 2 = 1 x 1 0 3 the 

expected number of events is 226. For = lOOGeV and |V),n | 2 = 1 x 10~ 3 the expected 
number of events is 4.4, while for = 500 GeV and | V ?i N| 2 = 1 it is 6.9. 

We see no evidence for a significant excess in data beyond the backgrounds predicted from the 
SM and we set 95% confidence level (CF) exclusion limits on the cross section times branching 
fraction for pp —> N jA~ —? jA-fi-qq as a function of m N , using the CL S method |49 tI 5T1 based on 
the event yields shown in Table [3] We use Poisson distributions for the signal and log-normal 
distributions for the nuisance parameters. The limits obtained are shown in Fig. [4] We do not 
consider the mass range below = 40 GeV because of the very low selection efficiency for 
the signal in this mass region. 

The behavior of the limit around = 80 GeV is caused by the fact that as the heavy Majorana 
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5 Results and discussion 


19.7 ft)- 1 (8 TeV) 



19.7 fb' 1 (8 TeV) 



19.7 fb' 1 (8 TeV) 



j_i ± |a=jj invariant mass (GeV) 


Figure 2: Kinematic distributions for the low-mass region after all selections are applied except 
for the final optimization requirements: leading muon pj (top), trailing muon pj (middle), 
and p-p-)\ invariant mass (bottom). The plots show the data, backgrounds, and two choices 
for the heavy Majorana neutrino signal: m\j = 40 GeV, |V)/n| 2 = 5 x 10 5 and = 80 GeV, 
|V fiN | 2 = 1 x 10-3. The backgrounds shown are from misidentified muons and from diboson 
(W), Higgs boson, triboson (WV), and ttW production. 
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Table 3: Observed event yields and estimated backgrounds for each Majorana neutrino mass 
point ()Hn). The background predictions from the irreducible SM backgrounds (SM Bkgd.), 
misidentified muon background (Misid. Bkgd.) and total background (Total Bkgd.) are shown 
in columns 2-4, while column 5 shows the number of events observed in data. The uncertain¬ 
ties shown are respectively the statistical and systematic components. As discussed in the text, 
the search was divided into low-mass (»zn < 80 GeV) and high-mass (zzz\j > 80 GeV) regions, 
and different selection criteria were used in the two regions. 


7/?n SM Bkgd. Misid. Bkgd. Total Bkgd. N 0 b s 

(GeV) 


40 

3.0 

± 

0.4 

± 

0.6 

6.7 

± 

0.9 

± 

1.9 

kO 

00 

± 

1.0 

± 

2.0 

7 

50 

3.0 

± 

0.4 

± 

0.6 

6.7 

± 

0.9 

± 

1.9 

9.8 

± 

1.0 

± 

2.0 

7 

60 

3.0 

± 

0.4 

± 

0.6 

6.7 

± 

0.9 

± 

1.9 

9.8 

± 

1.0 

± 

2.0 

7 

70 

3.0 

± 

0.4 

± 

0.6 

6.7 

± 

0.9 

± 

1.9 

9.8 

± 

1.0 

± 

2.0 

7 

80 

3.0 

± 

0.4 

± 

0.6 

6.7 

± 

0.9 

± 

1.9 

9.8 

± 

1.0 

± 

2.0 

7 

90 

8.7 

± 

0.7 

± 

1.7 

12.6 

• i 

: 1.1 

± 

: 3.5 

21.3 

.± 

: 1.3 

■ i 

: 3.9 

19 

100 

8.7 

± 

0.7 

± 

1.7 

11.7 

■± 

: 1.0 

i± 

: 3.3 

20.4 

: i 

: 1.2 

,± 

: 3.7 

19 

125 

7.9 

± 

0.6 

± 

1.5 

5.9 

± 

0.7 

± 

1.6 

13.8 

;± 

0.9 


: 2.2 

8 

150 

6.4 

± 

0.5 

± 

1.2 

3.6 

± 

0.6 

± 

1.0 

9.9 

± 

0.8 

± 

1.6 

7 

175 

4.4 

± 

0.4 

± 

0.8 

1.6 

± 

0.4 

± 

0.5 

6.0 

± 

0.6 

± 

1.0 

7 

200 

3.4 

± 

0.4 

± 

0.7 

0.8 

± 

0.3 

± 

0.2 

4.2 

± 

0.5 

± 

0.7 

5 

250 

1.9 

± 

0.3 

± 

0.3 

0.6 

± 

0.2 

± 

0.2 

2.5 

± 

0.3 

± 

0.4 

3 

300 

0.9 

± 

0.2 

± 

0.2 

0.1 

± 

0.2 

± 

0.0 

1.0 

± 

0.3 

± 

0.2 

1 

350 

0.9 

± 

0.2 

± 

0.2 

0.1 

± 

0.2 

± 

0.0 

1.0 

± 

0.3 

± 

0.2 

1 

400 

0.9 

± 

0.2 

± 

0.2 

0.1 

± 

0.2 

± 

0.0 

1.0 

± 

0.3 

± 

0.2 

1 

500 

0.9 

± 

0.2 

± 

0.2 

0.1 

± 

0.2 

± 

0.0 

1.0 

± 

0.3 

± 

0.2 

1 


neutrino gets close to the W-boson mass from below or above, the muon produced together 
with the N or the muon from the N decay have low pj, respectively. We also set limits on 
the square of the heavy Majorana neutrino mixing parameter times branching fraction of the 
N to a W boson and a muon, |V),n| 2 B(N —* W~p ), as a function of ni\. Figure [5] shows the 
resulting upper limits on | V ?i N| 2 f3(N —> 'W ± p - ) as a function of m\j. Assuming the theoretical 
prediction for the branching fraction for N —> W~p , we can extract limits on | 2 - These 
limits are shown in Fig. [6] 


6 Summary 

A search for heavy Majorana neutrinos in p~p~}) events has been performed in proton-proton 
collisions at a center-of-mass energy of 8 TeV, using a data set corresponding to an integrated 
luminosity of 19.7 fb -1 . No excess of events beyond the standard model background prediction 
is found. Upper limits at 95% CL are set on | V fl ^\ 2 , as a function of heavy Majorana neutrino 
mass 77 Zn, where 1 /,,n is the mixing element of the heavy neutrino N with the standard model 
muon neutrino. 

These are the first direct upper limits on the heavy Majorana neutrino mixing for > 200 GeV. 
The limits for m n between 45 GeV and 90 GeV are comparable to the best limits from LEP Il20l 
[2T| , which are derived from Z —> 17 N and are restricted to masses below approximately 
90 GeV. The limits reported here extend well beyond this mass and are significantly better 
than previous LHC limits obtained at a center-of-mass energy of 7 TeV. For = 90 GeV 
we find \V^\ 2 < 0.00470. Furthermore, the direct limits on \%\jj have been extended from 
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6 Summary 


OTn « 200 GeV up to « 500 GeV. At = 200 GeV the limit is IV^nI 2 < 0.0123, and at 
Wn = 500 GeV the limit is |V;<n| 2 < 0.583. 
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19.7 fb 1 (8 TeV) 



19.7 fb' 1 (8 TeV) 



19.7 fb' 1 (8 TeV) 



Figure 3: Kinematic distributions for the high-mass region after all selections are applied except 
for the final optimization requirements: leading muon pj (top), trailing muon pj (middle), 
and p-p-)} invariant mass (bottom). The plots show the data, backgrounds, and two choices 
for the heavy Majorana neutrino signal: = 100 GeV, | | 2 = 0.001 and = 300 GeV, 

|V fiN | 2 = 0.05. The backgrounds shown are from misidentified muons and from diboson (VV), 
Higgs boson, triboson (VVV), and ttW production. 
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6 Summary 



m (GeV) 

N v ' 


Figure 4: Exclusion region at 95% CL in the cross section times branching fraction, as a function 
of the heavy Majorana neutrino mass. The long-dashed black curve is the expected upper limit, 
with one and two standard deviation bands shown in dark green and light yellow, respectively 
The solid black curve is the observed upper limit. The region above the exclusion curves is 
ruled out. 
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Figure 5: Exclusion region at 95% CL in the square of the heavy Majorana neutrino mixing 
parameter times branching fraction of the N to a W boson and a muon, as a function of the 
heavy Majorana neutrino mass: (| V ? ,n| 2 !3(N —> ) vs. m N ). The long-dashed black curve 

is the expected upper limit, with one and two standard deviation bands shown in dark green 
and light yellow, respectively The solid black curve is the observed upper limit. The regions 
above the exclusion curves are ruled out. The lower panel shows an expanded view of the 
region 40 GeV < < 250 GeV. 
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6 Summary 



50 100 150 200 250 300 350 400 450 500 

m (GeV) 

N V ' 



m N (GeV) 

Figure 6: Exclusion region at 95% CL in the square of the heavy Majorana neutrino mixing 
parameter as a function of the heavy Majorana neutrino mass: ((V^nI 2 vs. The long- 

dashed black curve is the expected upper limit, with one and two standard deviation bands 
shown in dark green and light yellow, respectively The solid black curve is the observed upper 
limit. Also shown are the upper limits from other direct searches: L3 |[20| , DELPHI f2lT| , and 
the upper limits from CMS obtained with the 2011 LHC data at yjs = 7TeV 11221 . The regions 
above the exclusion curves are ruled out. The lower panel shows an expanded view of the 
region 40 GeV < < 250 GeV. 
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